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ABSTRACT 

We present a kinematic analysis of the globular cluster (GC) system of the giant ellip- 
tical galaxy NGC 4365 and find several distinct kinematic substructures. This analysis 
is carried out using radial velocities for 269 GCs, obtained with the DEIMOS instru- 
ment on the Keck II telescope as part of the SAGES Legacy Unifying Globulars and 
Galaxies Survey (SLUGGS). We find that each of the three (formerly identified) GC 
colour subpopulations reveal distinct rotation properties. The rotation of the green 
GC subpopulation is consistent with the bulk of NGC 4365's stellar light, which 'rolls' 
about the photometric major axis. The blue and red GC subpopulations show 'nor- 
mal' rotation about the minor axis. We also find that the red GC subpopulation is 
rotationally dominated beyond 2.5 arcmin (^ 17 kpc) and that the root mean squared 
velocity of the green subpopulation declines sharply with radius suggesting a possible 
bias towards radial orbits relative to the other GC subpopulations. Additionally, we 
find a population of low velocity GCs that form a linear structure running from the 
SW to the NE across NGC 4365 which aligns with the recently reported stellar stream 
towards NGC 4342. These low velocity GCs have g' — i' colours consistent with the 
overall NGC 4365 GC system but have velocities consistent with the systemic ve- 
locity of NGC 4342. We discuss the possible formation scenarios for the three GC 
subpopulations as well as the possible origin of the low velocity GC population. 



1 INTRODUCTION 

Globular clusters (GCs), as some of the oldest and densest 
stellar systems in the universe, present intriguing questions 
about their own formation while also being very useful trac- 



ers of the formation of their host galaxy ( Brodie & Strader 



20061. They probably formed during violent star formation 



episodes in their host galaxy (or host galaxy progenitors) 
and therefore we can use their spatial distribution and kine- 
matic properties to unravel the formation history of galaxies. 
GCs are more numerous and more easily identifiable in el- 
liptical galaxies (where bright star forming clumps and dust 
obscuration is minimal) than in spiral galaxies. 

NGC 4365 is a giant elliptical galaxy (E3), Mb = -21.3 
Mv = —22.7 mag (Ferrarese et al. 2006 Kormendy et al. 



2009| , with various intriguing and unusual properties as 



highlighted below. We use these unusual properties of NGC 
4365 to constrain galaxy and GC formation scenarios by 



placing the extra information contained in this system in 
the overall framework of 'normal' giant elliptical galaxies. 

The GC systems of most giant elliptical galaxies are bi- 
modal in optical colours ( Brodie fc Strader|2006 Peng et al. 
|2006[ ). There is a growing body of evidence that the relative 
age difference between GC subpopulations is small ( Strader 



et al. 2005 1 and the distinct colours of the commonly ob- 



served subpopulations are due to bimodality in GC metal- 



licity ([Strader et al.||2007| [Woodley et al.|20"l0{ |Alves-Brito 



et al.|20lT I that are likely explained by two formation mech- 
anisms, sites or epochs. Possible formation scenarios for bi- 
modal metallicity distributions in GC systems are the major 



merger ( Zepf fc Ashman|1993 1, multiphase collapse (Forbes, 



Brodie & Grillmair 19971) and accretion scenarios ( Cote et al 



19981. An alternative view is that the bimodal GC colour 



distributions can be observed from a unimodal GC metal 
licity distribution, because of strong nonlinearity in the re- 
lationship between colour and metallicity ( Yoon et al.|2006 
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Figure 1. The positions of the six DEIMOS sUtmasks plotted on the Subaru Suprime-Cam (S-Cam) i' filter image. The scale is shown 
in the bottom left corner of the 35 X 27 arcmin image. It is centred on 0=12:24:26.824; 5=+07:19:03.52 (J2000.0). At a distance of 
23.1 ± 0.8 Mpc jBlakeslee et al.|2009[ l 1 arcmin = 6.72 kpc. 



Cantiello fc Blakeslee'2007"Blakcslce et al.|2010l ) and in this 
case the kinematic properties of each colour subpopulation 
are not expected to be significantly different. 

The GC system of NGC 4365 has three GC subpopu- 
lations, that includes an additional subpopulation of 'green' 
GCs between the commonly observed blue and red GC sub- 
populations. This additional subpopulation, thought at the 
time to be younger than the 'normal' two subpopulations, 



was discovered by Puzia et al. ( 2002 \ using a combination 
of near-IR and optical photometry. [Brodie et alT] ( |2005[ ) sug- 
gested that the three subpopulations might be found as tri- 
modality in optical colours and Larsen et al. (12005) found 



that the green subpopulation was restricted to small galac- 
tocentric radii. 

Specifically for NGC 4365, two independent methods 
have confirmed the GCs for all three subpopulations (blue, 
green and red) to be older than ~ 10 Gyr. Brodie et al. 



1 2005 1 observed a spectroscopic sample of 22 GCs to mea- 



sure their ages and metallicities using Lick indices from the 
Low-Resolution Imaging Spectrograph (LRIS) on the Keck 
telescope. They found that NGC 4365 has three old GC 
subpopulations with metal poor (blue), intermediate metal- 



licity (green) and metal rich (red) GCs. Chies-Santos et al. 
(2011 1 presented a near-IR photometry and optical analysis 



sible, and determined that there is no observable offset in 
the mean age between the GCs of NGC 4365 and GC pop- 
ulations of other giant ellipticals. This indicates that there 
is no significant population of young GCs in NGC 4365. 

In addition to its almost unique GC system, NGC 4365 
has very unusual stellar kinematic properties. The kinemat- 
ically distinct core (KDC) at its centre is relatively uncom- 
mon, seen in ~ 10 per cent of early type (E and SO) galaxies 
in the ATLAS'^'' volume limited survey (Krajnovic et al. 
|2011[ ). Also very uncommon is the starlight kinematics out- 
side the KDC which 'rolls' about the major axis rather than 
rotating about the minor axis ( Surma &: Bender|1995 1. Less 
than 5 per cent of galaxies, in the ATLAS'^^ volume limited 



survey, show this kinematic behaviour. Davies et al. (20011 



used the SAURON integral field spectrograph to map its 
kinematic and metallicity structure in two dimensions, find- 
ing stars both inside and outside the KDC to be older than 
10 Gyr. NGC 4365 is one of only two early type galaxies, in 
the ATLAS^-° sample of 260, that have both a KDC and a 
near 90'' misalignment between the photometric and kine- 
matic major axes (i.e. rolling rather than rotating stars). 



Lastly Bogdan et al. (20121, see also Mihos et al. (2012, 



of 99 GCs with much higher precision than previously pos- 



in prep), presented evidence of a stream of stellar light ex- 
tending ~ 200 kpc southwest and ~ 100 kpc northeast from 
NGC 4365. This strong indicator of a very recent ~ 10 : 1 
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merger (occurring with in the last few Gyr) presents a pos- 
sible cause for the unique properties of NGC 4365. 



[Fe/H] 



Blom et al. (20121 presented an optical photometry 



analysis of ~ 4000 NGC 4365 GCs, using the most spa- 
tially extended sample of GCs analysed to date (observed 
with Subaru/Suprime-Cam and supplemented with archival 
imaging from Hubble Space Telescope/ Advanced Camera 
for Surveys) and found different spatial distributions, sizes 
and mass distributions for the three subpopulations of NGC 
4365 GC system. They also concluded that separate forma- 
tion scenarios are required to explain the existence of the 
three separate subpopulations of GC in NGC 4365. 

Kinematic analysis of the three GC subpopulations in 
NGC 4365 is key to disentangling the possible formation 
scenarios of GC systems. Several recent investigations have 
shown that the kinematic features of the two standard GC 
subpopulations (seen in other giant elliptical galaxies) are 



different (Lee et al. 2010 Arnold et al. 2011 Foster et al 



2011). When combined with galaxy formation models this 



information constrains the possible formation scenarios of 



GCs. For example, Foster et al. (20111 compare kinematics 



of the GC subpopulations in NGC 4494 with simulations 



such as those described in Bekki et al. ( 2005 1 and conclude 



that the galaxy has undergone a recent major merger of sim- 
ilar disk galaxies. The addition of the kinematic properties 
of a third GC subpopulation limits the possible formation 
scenarios even further. 

In Section 2 we first describe the preparation, observa- 
tion and reduction of the spectroscopic sample of GCs and 
then investigate the colour/metallicity and line-of-sight ve- 
locity distributions of the GCs. The separation of the GCs 
into subpopulations and analysis of the kinematic features 
of the three subpopulations are presented in Section 3. We 
then discuss the results and summarise our conclusions re- 
garding the possible formation scenarios for NGC 4365 in 
Sections 4 and 5. The Appendices contain the kinematic fits 
for all the different methods of subpopulation separation we 
investigated, and an alternative assumption of position an- 
gle for kinematic fitting. 



2 SPECTROSCOPIC SAMPLE 
2.1 Data acquisition 

The GC spectroscopic sample presented in this work was 
obtained from six multi-object DEIMOS slitmasks observed 
in clear dark time on the Keck II telescope. This galaxy 
was observed as part of the SAGES Legacy Unifying Globu- 
lars and Galaxies Survey (SLUGG^ (Brodie et al. 2012 in 
prep). Four masks were observed on 2010, January 11 and 
12 with seeing between 0.6 and 0.9 arcsec and a further two 
masks were observed on 2010, February 18 with seeing of 
~ 0.85" 



Blom et al. (20121 used Subaru/Suprime-Cam (S- 



Cam) and Hubble Space Telescope / Advanced Camera for 
Surveys (HST/ACS) imaging to identify GC candidates as- 
sociated with NGC 4365 and publish extinction corrected g' , 
r' and i' magnitudes for S-Cam and/or g and z magnitudes 
for HST/ACS candidates. A total of 443 GC candidates with 
i' brighter than 23 mag were placed in slits, for which we 

http://sluggs.swin.cdu.au/ 
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Figure 2. Colour distribution (reddening corrected) of GCs with 
radial velocity measurements compared with the total GC photo- 
metric sample. The grey histogram shows the kinematic sample 
and the white histogram shows the distribution of photometric 
candidate GCs brighter than the turnover magnitude, i' = 23.6, 
and over all observable radii, ~ 20 arcmin, (scaled down by a 
factor of 2.7 for comparison purposes). The kinematic sample is 
a subsample of the photometric candidate sample. The top x- 
axis shows the [Fc/H] metallicity determined from the empirical 
conversion published in |Lee et al.| l [2010[ |. We note the three GC 
subpopulations with peaks at g' — i' = 0.8, 0.97 and 1.13. We 
also note that the blue GCs in the photometric candidate sample 
are underrepresented in the kinematic sample with respect to the 
number of redder GCs {g' — i' > 0.9). 



aimed to obtain a spectrum of high enough signal-to-noise 
(S/N ~ 5) to determine a line-of-sight radial velocity. For 
all six masks a 1200 1 mm~^ grating was used, with 1 arcsec 
wide slits and centred on 7800 A to obtain ~ 1.5 A resolu- 
tion. This set up (see also Foster et al.||2011 Arnold et al 



201 1| [Strader e t al. 2011, Romanows ky et al.||2009| |2012 l 
allows observations from ~ 6550 - 8900 A , which covers the 
wavelength range in which we expect the red-shifted Cal- 
cium II triplet (CaT) absorption features (8498, 8542 and 
8662 A). It sometimes also includes the red-shifted Ha line 
(6563 A) at the blue end of the observed spectrum. 

Footprints of the six DEIMOS masks plotted on the 
S-Cam i' filter image, used to identify GCs around NGC 
4365, are shown in Fig.[T] We obtained total exposure times 
of 110, 40, 105, 96, 120 and 75 minutes for each of the six 
masks respectively. For each mask the total exposure time 
was split into three or four separate exposures to minimize 
the effects of cosmic rays. The total exposure times for each 
mask were independently constrained by seeing and object 
visibility during the night. The median seeing was 0.8, 0.7, 
0.75, 0.7, 0.9 and 0.7 arcsec for each of the six masks respec- 
tively. Masks two and six were observed at the end of their 
respective observing nights. 

The raw data were reduced to one dimensional spec- 
tra using a customised version of the DEEP2 galaxy survey 
data reduction pipeline (IDL SPEC2D; [Cooper et al.|2012 
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Newman et al.|["2012 |. The pipehne uses dome flats, NeAr- 
KrXe arc lamp spectra and sky hght visible in each slit to 
perform flat fielding, wavelength calibration and local sky 
subtraction, respectively. It outputs the object spectrum as 
well as the locally subtracted sky spectrum. 



2.2 Obtaining line-of-sight velocities 

We obtained reliable line-of-sight radial velocities for 252 
GCs. This increases the previously available NGC 4365 spec- 



troscopic GC sample (33 individual GCs in Larsen et al, 
|2003[ [Brodie et al.|[2005{ combined) by an order of magni- 
tude. 

The line-of-sight radial velocity is determined by the 
peak of the cross correlation function (obtained using the 
IRAF task RV.FXCOR) between the spectrum and 13 stel- 
lar templates. The velocity measurement of a GC is consid- 
ered reliable if the correlation function has a single, identifi- 
able peak and at least two of the four absorption lines (three 
CaT and one Ha) are clearly visible in the red-shifted spec- 
trum. In a few cases one or other of these visual checks was 
not convincing and the resulting GC velocity is recorded as 
a borderline velocity measurement (17 objects in total). The 
final catalogue of 269 radial velocities for NGC 4365 GCs are 
given in Pota et al. (2012, in prep.) 

The colour distribution of the kinematic spectroscopic 
sample is compared with that of the photometric sample 
in Fig. |2] The photometric sample includes GC candidates 
with g' — i' colours between 0.5 and 1.4, showing overdensi- 
ties at g' — i' 0.8 and between 0.95 and 1.2 as well as an 
underdensity at g' — i! ~ 0.9. Ideally the kinematic sample 
will have the same colour distribution, scaled down in num- 
ber, as the photometric sample. In this ideal case we would 
be able to extrapolate conclusions drawn from the spectro- 
scopic samples to the entire GC population without caveat. 
Fig. [2] shows that while we do not completely achieve this 
ideal, the spectroscopic sample does cover the full colour 
range and shows a significant underdensity at intermedi- 
ate colours between the blue and the green subpopulations 
(seen at g' — i! ~ 0.9). In the spectroscopic sample blue GCs 
are under represented with respect to the green/red GCs 
compared to the photometric sample. This is because the 
six DEIMOS pointings have a much smaller radial extent 
(~ 12 arcmin) than the S-Cam imaging (~ 20 arcmin) and 
at larger galactocentric radii there are many more blue GCs 
than red or green GCs (see Fig, [s] and [Blom et al.|20"T2l ). 

We expect to see a roughly Gaussian distribution in the 
line-of-sight velocities of GCs around NGC 4365. In Fig. [3] 
we show a Gaussian distribution overplotted on the velocity 
distribution of NGC 4365's GCs. The peak is set to the 
literature val ue for NGC 4365' s recession velocity [vsys ~ 
1243 km s~ MSmith et al.|2000 l and the standard deviation 
is set to 237 km s~ . This standard deviation was calculated 
from one side of the velocity distribution, i.e. from 68 percent 
of the GCs with Vioa > 1243 km s~^. The peak agrees well 
with the observed GC velocity distribution but an excess of 
GCs is seen above the overplotted Gaussian at low velocities 
in Fig.m 

We can calculate the skewness (^3) of the velocity sam- 
ple using the standardised third moment (Joanes & Gill 
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Figure 3. Linc-of-sight velocity distribution for GCs in the kine- 
matic sample. The distribution of the secure line-of-sight velocity 
measurements (252 GCs) is plotted with a solid line and the bor- 
derline velocity measurements (17 GCs) are added to the distribu- 
tion with a dashed line. The plotted Gaussian distribution has a 
peak at 1243 km s~^ and a velocity dispersion of 237 km s~^. The 
velocity distribution of the GCs agrees well with the Gaussian at 
high velocities but is skewed to low line-of-sight velocities. 



19981, 



6 



Vjv(jv-i) Eti(^»-^) 

N{N - 2) s3 



(1) 



where is the total number of GCs in the sample, Xi is the 
line-of-sight velocity of the z*'' GC, x is the mean line-of- 
sight velocity of the sample and s is the standard deviation 
of the sample. The calculated skewness is not strongly de- 
pendent on whether the entire GCs sample (including bor- 
derline velocity measurements) or only the secure velocity 
measurements are used, 

6,aii = -0.20 

^3, secure — 0.19 

The standard error in skewness {e^^ ) can be calculated from 
the number of GCs in the sample ( Joanes fc Gill||1998^ , 



«C3 = 



6A(A- 1) 



(A-2)(A + l)(A + 3) 



= 0.15 



(2) 



The measured negative skewness of the GC line-of-sight ve- 
locity distribution (—0.19 ± 0.15) is only marginally signifi- 
cant, given the relatively small number of GCs contributing 
to the skewness. We further investigate the possible causes 
of this negative skewness in the velocity distribution. 



2.3 Low velocity GCs 

We explore the possibility that the GCs causing marginal 
skewness in the measured velocity distribution are associ- 
ated with the stream of stars that were found crossing NGC 
4365 (North East to South West) from analysis of wide-field, 
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Figure 5. GC surface density with colour and galactocentric radius. The colour scale indicates object number per arcmin'^ per colour bin 
ranging from black at low density to red at high density. The colour-radius position of spectroscopically confirmed GCs are overplotted with 
white dots. GCs with radial velocities more than 2a from the host galaxy are marked with black plusses (positive) and dots (negative). 
The distribution of GCs with observed line-of-sight radial velocities is consistent with the underlying distribution of photometrically 
observed GCs in the same radial range and there is no indication that the GCs with > 2cf low velocities deviate from the underlying 
distribution either. 
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Figure 4. Line-of-sight radial velocity {Vioa) of GCs plotted 
against galactocentric radius. Filled black circles show the 252 
GCs with reliable radial velocity measurements, while grey cir- 
cles show borderline GC detections. The solid black line is plot- 
ted at the galaxy systemic velocity Vsys = 1243 km s~^, dotted 
lines show la and 2cr envelopes and the dashed line shows the 3(t 
envelope. The a values were calculated in bins of 29 GCs with 
Vlos > 1243 km s^^, excluding borderline GCs. Note the large 
number of low velocity {Vios ~ 700 km s^^) > 2a outliers. The 
solid grey line is plotted at the syste mic velocity of the nearby 
galaxy NGC 4342 {Vsya = 751 km s~ MGrogin et al.|l998| . which 
lies ~ 35 arcmin from NGC 4365 



very deep B filter imaging fBog dan et al.|2012 Milios et al. 
2012, in prep). The stream seems to be an indication of an 
ongoing merger with a small galaxy. The distinct velocities 



of these GCs may indicate that they are not associated NGC 
4365's GC system and the recently discovered stream is a 
likely candidate. Fig. [1] shows the GC line-of-sight veloci- 
ties plotted against galactocentric radius. To define which 
GC velocities are outliers we used the GC velocities > 1243 
km s~^ (where the velocity distribution is not affected by 
skewness) to calculate the standard deviation of the sample, 
a. The calculation of cr was done in radial bins of 29 GCs 
each (excluding borderline velocity measurements) . The cal- 
culated la, 2a and 3a envelopes are plotted above and below 
the systemic velocity and show clearly that there is a larger 
number of GCs with low velocities outside the 2cr envelope 
(13) than with high velocities (3). 

When we plot these velocity outliers on a map of GC 
surface density with colour and galactocentric radius there 
is no strong indication that the low velocity outliers follow 
a separate distribution to the colour-radius distribution of 
NGC 4365 GCs (see Fig. |5|. The density map shows the 
distribution of photometrically identified GCs overplotted 
with the kinematic sample of GCs (those with reliable ra- 
dial velocities). The blue and red subpopulations, as well 
as an overdensity of green GCs in the central regions, are 
visible in the density plot and the GCs with radial veloc- 
ity measurements follow the colour-radius distribution well. 
The low velocity outliers tend to have red or green colours 
in the inner regions and blue colours in the outer regions of 
the galaxy, which is consistent with the general distribution 
of NGC 4365 GCs in colour-radius space. 

In Fig. |6] the positions of the low velocity GCs are plot- 
ted on a map of the spatial distribution of the GCs with 
line-of-sight velocities. All but two of the low velocity GCs 
form a line running NE to SW across NGC 4365. This spatial 
structure of low velocity GCs is aligned with the recently dis- 
covered tidal stream ( Bogdan et al.|2012 Mihos et al. 2012, 
in prep) and would seem to indicate that the majority of the 
low velocity GCs are associated with the stream. 

In summary, the evidence does not conclusively show 
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Figure 6. Spatial distribution of GCs showing individual radial 
velocities. The colour of individual points correspond to their ra- 
dial velocity and the scale runs from 500 (dark blue) to 1800 km 
s~^ (dark red). GCs writh radial velocities more than 2a from the 
host galaxy are marked with plusses (positive) and dots (nega- 
tive) . 



that the low velocity GCs are associated with the stellar 
stream. Their distribution in colour-radius space is not dis- 
tinguishable from the overall distribution of NGC 4365 GCs 
but their spatial distribution does suggest that they are dif- 
ferent from the bulk of NGC 4365 GCs. We exclude all GCs 
with line-of-sight velocities that lie outside the 2a" envelopes 
and GCs with borderline velocity measurements from fur- 
ther kinematic analysis. This conservative sample ensures 
that the kinematic results are unaffected by unreliable GC 
velocity measurements and that the sample is not skewed 
to low velocities. We note here that including the low ve- 
locity outliers makes no significant difference to subsequent 
kinematic results, except to increase the uncertainties in our 
fits. 



3 KINEMATICS OF THE GC 
SUBPOPULATIONS 

3.1 Kinematic model description 

To investigate the detailed kinematic properties of the three 
GC subpopulations we fit an inclined disk model to the GC 
system of NGC 4365, using the individual GCs as tracers. 
This inclined disk model is able to fit for rotation velocity 
and the position angle of that rotation with our sparse and 
non-uniform sampling of tracer particles. To fit this model 
we minimize (as in Foster et al.|20lT |; 



{Vob 



0-2 + (AK, 



-I- ln(cr + {/Wobs,i) , 



Vrr, 



(3) 



(4) 



1 + 



tan(PAi-PAfci„) 
Ikin 



where Vobs,i is the observed line-of-sight velocity, 
AVobs.i is the uncertainty on the observed velocity, a is the 
fitted velocity dispersion of the sample, Vmod,i is the mod- 
elled velocity given by equation (4), Vaya is the recession 
velocity of NGC 4365, Vrot is the rotation velocity fitted to 
the sample, PAi is the position angle of the i^^ GC, PAfei„ is 
the kinematic position angle fitted to the sample and Qkin is 
the kinematic axis ratio (b/a) fixed for the sample. In con- 
trast to Foster et al. (20111 we do not use a minimisation 



algorithm to determine the minimum value of A but eval- 
uate the function for all reasonable values of Vrot, o and 
PAfcin. This brute force method ensures that we have found 
the global minimum for the function and therefore the best 
fitting values. 

We use the bootstrap Monte Carlo method to obtain 
an estimate of the uncertainty of the model fit. A ran- 
dom sample of GCs is chosen from the original GC sam- 
ple (equal in number to the original sample but allowing 
the data from a given GC to be used more than once) and 
then A is minimized for that randomised sample. This is 
repeated 1000 times and the uncertainties are determined 
by calculating the central 68 per cent confidence intervals 
of the distribution (la) for each of Vrot, a and PAfci„. To 
use a brute force minimisation method for each iteration 
would be prohibitively time consuming and therefore we use 
a PoweU minimisation algorithm (from SCIPY. OPTIMIZE 
in the PYTHON framework) with initial estimates set to 
the solution obtained from the kinematic best fitting values. 
The Powell minimisation algorithm is marginally sensitive 
to the initial conditions but we found no significant differ- 
ence in the uncertainties obtained using this method and 
those found for a test case where the brute force method 
was used instead. 

In general the la uncertainties are not symmetric about 
the best fitting value because the distribution obtained from 
the bootstrap algorithm is not symmetric about the peak 
value. We calculate a central 68 per cent confidence interval 
rather than a symmetric one. The peak value of the distri- 
bution aligns well with the best fitting value in most cases. 
We do not artificially align the peak value of the distribu- 
tion with the best-fit value and consequently it is possible 
for the best fitting value to fall outside the la range of the 
distribution. This is a clear indication that for such a case 
the model is inappropriate for the data. 



3.2 Kinematics as a function of colour 

We use our kinematic model to analyse the kinematics of 
NGC 4365's GCs as a function of colour. The GCs with 
velocities less than 2a away from the galaxy systemic ve- 
locity and secure velocity measurements are used to form a 
rolling sample of 30 GCs as follows. The 30 GCs with bluest 
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Figure 7. Kinematics for NGC 4365 GCs as a function of g' — i' 
colour. Kinematics are fitted for a rolling sample of 30 GCs and 
are plotted at the mean colour of the sample. The top panel shows 
the fitted rotation velocity {Vrot), the middle panel shows the 
fitted velocity dispersion {a) and the bottom shows the fitted 
kinematic position angle (PAj.i„). Solid lines show the kinematic 
fit to the data and the grey filled area shows the la uncertainties 
(obtained from 1000 bootstrap trials). On the bottom panel the 
photometric major axis of NGC 4365 (Blom ct al. 2012| is plotted 
with dashed lines and the minor axis is plotted with dotted lines. 
The bottom panel shows that the very blue, green and very red 
GCs all have constrained position angles of rotation and that the 
blue and red subpopulations rotate in a different direction to the 
green subpopulation. Where there is significant contamination 
between subpopulations in colour the rotation is unconstrained. 



colours form the first subsample. To form the next subsam- 
ple the bluest GC is discarded while the GC with the bluest 
colour that was not included in the first subsample (i.e. the 
next reddest) is added and this process is repeated until all 
GCs have been included in a subsample. We chose the sam- 
ple size of 30 as a best compromise between samples large 
enough to robustly constrain the kinematics and samples 
small enough to give adequate colour resolution. The kine- 
matic fitting procedure is carried out on each subsample and 
the results are shown in Fig. [7] In this case we fix the kine- 
matic axis ratio to one {qkin ~ 1) as the photometric axis 
ratio for the galaxy light as well as blue, green and red GC 
subpopulations are all different {qbiue = 0.56, Qgreen = 0.70, 
Qred ~ 0.97 and qgai — 0.75 from Blom et al. 



2012| . This 



assumption of qkin = 1 is equivalent to assuming a circular 
rotation disk and is chosen as the simplest case scenario. 

Fig. [7] shows that the GCs of NGC 4365 have signif- 
icantly different kinematic features, which depend on GC 
colour. The changes in position angle seen in Fig. [7] show 
this very clearly, e.g. blue GCs rotate with PAkin that is 
consistent with the photometric major axis, green GCs ro- 
tate with PAkin consistent with the photometric minor axis 
and red GCs with rotate with PAkin that is between the pho- 
tometric major and minor axes of NGC 4365. At the colours 
where we expect significant overlap between subpopulations 
the rotation cannot be constrained. Here the values fitted 



for the PAkin become unstable while the values fitted for 
Vrot drop to zero and become inconsistent with the errors 
calculated. Where it is not possible to constrain the posi- 
tion angle of the rotation, the amplitude of rotation is also 
unconstrained. 



Foster et al. (20111 noted that Vrot is artificially in- 



creased when the PAkin is allowed to vary. jStrader et al.| 
(2011 1 describe this artificial increase as a rotation bias and 



note that this bias is important when the amplitude of ro- 
tation is low, e.g. 



a V jV 



(5) 



Here A'^ denotes the number of GCs in the sample. When we 
apply this criterion to the kinematic fit to NGC 4365's GCs 
as a function of colour the colour ranges where the PAkin 
is weU constrained {g' - j' ~ 0.79 - 0.84, 0.97 - 1.02 and 
1.12 — 1.17) the rotation parameter {Vrot/o') is well above 
the criterion value for this sample size (A'' — 30). We also 
see that the rotation parameter fails the criterion (i.e. is 
affected by rotation bias) in areas where the PAkin is poorly 
constrained (e.g. 0.87 < g' — i' < 0.95) even when the fitted 
Vrot is measured to be significantly greater than zero. 

In the blue, green and red colour ranges where we do 
constrain the kinematics the blue subpopulation Vrot and a 
ranges are 50 — 140 and 140 — 160 km s~^ respectively while 
the green subpopulation Vrot and a ranges are 80 — 200 
and 180 — 220 km s~^ respectively. The rotation of the red 
subpopulation is range is Vrot ~ 80 - 180 km s"^ while 
its velocity dispersion decreases with increasingly red GC 
colour, a ~ 220 to ~ 200 km s"\ 

This observation of distinct kinematic features is a con- 
firmation that the three GC subpopulations of NGC 4365 



seen in colour and other photometric properties ( Larsen 
et al.|2003| [Brodie et al.|2005| |Blom et al.|2012[ ) are indeed 
separate. They are likely to have formed via different mecha- 
nisms or at different epochs. To better constrain the position 
angle (PAkin) and derive clearer values for the amplitude of 
rotation (Vrot) and dispersion (a), we need to divide the 
GCs into three uncontaminated representative samples of 
the three GC subpopulations. 



3.3 Dividing the sample into three subpopulations 



Blom et al. (20121 divided NGC 4365's GCs into three 



subpopulations according to a probability of belonging to 
that subpopulation. They assigned each GC a probabil- 
ity of belonging to the blue, green and red subpopulations 
[Pbiub + Pcreen + PRed ~ 1) based ou its colour and galac- 
tocentric radius. "The process used to determine the blue- 
green-red probability of a GC starts by determining the rel- 
ative number of GCs in each subpopulation at the galac- 
tocentric radius of the GC. This is done by comparing the 
values of the radial surface density profile for each subpopu- 
lation. The relative numbers of GCs in each population are 
used to scale the normalisations of the three Gaussian distri- 
butions and comparing the relative values of all three Gaus- 
sian distributions at the colour of the GC [they] calculate 
the probability of an object being blue, green, or red." ( Blom 
et al.|2012| ) Their high probability subsamples (P > 0.95) do 
not include any green GCs and therefore we define medium 
probability samples with P > 0.80 to represent the three 
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Table 1. Summary of kinematic fits to various divisions of the three GC subpopulations. We tabulate the rotation velocity {Vrot), 
velocity dispersion (ct) and kinematic position angle (PAj.i„) for each subpopulation in the medium probability (probability denoted 
with Pcolour) Jsnid best colour cut splits. Numbers of GCs in each subpopulation representative sample are given in brackets after the 
definition parameters. The major axis of NGC 4365's stellar light is 42° or 222°. 



Split 
Definition 


Vrot O" PAfci„ 

(kms-1) (kms-1) (°) 


Vrot 

(km s"-"^) 


(kms-1) (°) 


Vrot cr 

(km s-1) (km s-1) 


(°) 


medium probability 


PBlue > 0.8 (53) 
49±f| 179±i? 254±|| 


POreen > 0.8 (53) 

83±i 229±ii 145±i? 


Pfled > 0.8 (65) 
39±ff 218±it 


i74±ii 


colour cut 


g' -i' < 0.85 (61) 
67±ff 178±J|! 269±|| 


0.9 < g' 

81±i! 


- i' < 1.05 (77) 
217±10 144±i^ 


g' -i' > 1.1 (56) 
94±|| 202±i5 


194±|| 




Figure 8. Observed line-of-sight velocity is plotted against position angle for the blue (left), green (centre) and red (right) GC subpop- 
ulations. We show the subpopulations as defined by the best colour cut split (see Table [ij. The fitted rotation velocity is plotted with 
a solid sinusoid, and an envelope of twice the fitted velocity dispersion is shown with dotted lines for each subpopulation. The position 
angle (PAfci„) of rotation is 269°, 144° and 194° for the blue, green and red subpopulations respectively. The photometric major axis 
of the galaxy is 42/222° and the overdensity of blue GCs at ~ 210° is due to the strong elongation of the blue GC subpopulation along 
the major axis of the galaxy. 



subpopulations. Of GCs with radial velocities, 53 blue, 53 
green and 65 red GCs form part of the medium probability 
samples. 

The kinematic fits to the medium probability fits are 
shown in Table [l] The blue and green sample fits agree 
with the kinematic properties found in Fig. [7] but the fit 
to the red sample is not well constrained {PAkin varies by 
115°). The blue medium probability sample contains GCs 
with colours bluer than g — i! ^ 0.85 and the green sam- 
ple contains GCs with colours between g' ~ i' ~ 0.85 and 
g' ~ 1.05. These samples contain relatively few GCs but 
do not overlap and are mostly free of contamination from 
the other subpopulation. The red medium probability sam- 
ple contains more GCs but also contains GCs with colours as 
'blue' as — i' ~ 1.0 and is likely affected by contamination 
from green GCs. We expect some contamination particularly 
in the red sample because the medium probability red GCs 
extend almost as far as the peak colour of the green subpop- 
ulation ( 



i' = 0.1 



Blom et al. 20121. The probability 



that these GCs belong to the green subpopulation is close 
to 20 percent and not trivial. 

We also experimented with splitting the GCs into three 
subsamples at different g' — i' colours, detailed in Appendix 
A. To determine the most accurate kinematics for the three 
colour subpopulations of NGC 4365 we need to include as 
many GCs as possible in each sample while minimising the 



number of GCs from one subpopulation contaminating the 
sample for another subpopulation. The best kinematic re- 
sults are determined by the fit which has the minimum 
value of A/nd/ (see equation 3) and the most constrained 
PAfci„. Salient results are shown in Table [T] and details dis- 
cussed in Appendix A. The best kinematic fits are found for 
the 61 blue GCs with g' ~ i' < 0.85, 77 green GCs with 
0.9 < g' - i' < 1.05 and 56 red GCs with g' - i' > 1.1. 
The results for the medium probability and best colour cut 
samples are consistent within the uncertainties, but the best 
colour cut sample contains larger numbers of blue and green 
GCs and a relatively uncontaminated sample of red GCs. 
The red sample in this case ends close to the expected peak 
in the red subpopulation {g' — i' — 1.13). Henceforth we use 
the best colour cut sample for the subpopulation represen- 
tative samples. 

The kinematics are fitted twice for each colour split def- 
inition, once with the axis ratio fixed to the photometric 
value for each GC subpopulation {qkin ~ Qphotom) and once 
with the axis ratio fixed to unity {qkin = 1) as before. We 
do not have enough GCs in each sample to fit for qtin and it 
is not clear that the kinematic axis ratio would necessarily 
be equal to the photometric axis ratio, especially since Fig. 
[T] leads us to expect rotation in different directions to the 
photometric position angle. We do not find any significant 
differences in the fitted kinematics between the two axis ra- 
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tio cases and conclude that this kinematic fitting method is 
not very sensitive to the kinematic axis ratio (see Appendix 
A). All further kinematic fits are done in the case where 
Qkin = 1 as this is the more general case. 

We compare the rotation of these subpopulation repre- 
sentative samples (henceforth referred to as the blue, green 
and red subpopulations) with each other and with the kine- 
matics of the central starlight (to ~ 0.5 arcmin) as mea- 
sured by Davies et al. (20011 and Krajnovic et al. (20111. 





300 F 



galactocentric radius (kpc) 
20 30 



They found that the KDC rotates with a PAfei„ of 38 ± 5° 
aligned with the photometric major axis of the galaxy (42°) 
and that the bulk of the starlight rotates with a PAfei„ of 
145±6.5°, very close to the minor axis of the galaxy (132°). 
Using the Stellar Kinematics from Mu l tiple Slits (SKiMS ) 
technique (see also Proctor et al.|[2009 Foster et al.||20lT |, 
Arnold et al. (2012, in prep) find that the minor axis ro- 
tation of the starlight continues out to ~ 4.5 arcmin. We 
find that the blue subpopulation rotates with a PAfei„ be- 
tween 235° and 297°, not consistent with either the major 
(222°) or the minor (312°) axes of NGC 4365. The green 
subpopulation rotates with a PAfci„ of 144±23°, consistent 
with the minor axis of the galaxy and the rotation of the 
bulk of the starlight. The red subpopulation rotates with a 
PAfein between 170° and 217°, almost consistent with the 
major axis but in the opposite direction to the stars of the 
KDC. The kinematic properties that we have found by using 
the best representative samples for each GC subpopulation 
agree with the kinematics seen in Fig. [7] where GCs kine- 
matics are plotted as a function of colour. 

Fig. ?? shows the combined rotation velocity and ve- 
locity dispersion. These are plotted separately for each sub- 
population and the best-fit models show good agreement 
with the observed radial velocities of the subpopulations. 
The non-uniform azimuthal distribution of GCs, seen par- 
ticularly in the blue subpopulation, but also in the green 
subpopulation, is partly due to the elliptical spatial distri- 
bution of these subpopulations (axis ratio quue = 0.56 and 
Qgreen = 0.70, scc Blom ct al. 2012 1 and partly due to the 



asymmetric spatial coverage of the spectral observations. We 
do not expect this azimuthal bias to significantly influence 
the kinematic flts. 




2 3 4 5 

galactocentric radius (arcmin) 



Figure 9. Kinematics as a function of radius for NGC 4365's 
three GC subpopulations. GC subpopulations are defined using 
the best colour cut in Table [l] Kinematics are fitted for a rolling 
sample of 22 GCs and are plotted at the mean galactocentric 
radius of the sample. The top panel shows the fitted rotation 
velocity (Vrot), the middle panel shows the fitted velocity dis- 
persion (cr) and the bottom panel shows the fitted kinematics 
position angle (PAf;i„). The blue subpopulation (all GCs bluer 
than g' — i' = 0.85) is plotted with blue dashed lines, the green 
subpopulation (GCs between g' — i' = 0.9 and 1.05 with galacto- 
centric radius < 5 arcmin) is plotted with green dash-dotted lines 
and the red subpopulation (all GCs redder than g' — i' = 1.1) 
is plotted with red dotted lines. The lines show the kinematic fit 
to the data and the filled areas show the la uncertainties (ob- 
tained from 1000 bootstrap trials). The black points show stellar 
kinematics extracted along the minor axis of NGC 4365 from the 
SAURON data ( jde Zeeuw et al.||2002"l l. The effective radius of 
the galaxy starlight is plotted with a vertical line at 2.1 arcmin 
( |Blom et al.|2012] |. On the bottom panel the photometric major 
axis of NGC 4365 is plotted with dashed lines and the minor axis 
is plotted with dotted lines. The variation of the PApf^at is less 
than 3° and not visible on this scale. The kinematic fits to the 
three subpopulations show clearly different features. 



3.4 Radial kinematics for three subpopulations 

We use the best colour cut samples to analyse the variation 
of kinematics with galactocentric radius for each subpopu- 
lation. The samples contain 61 and 56 GCs for the blue and 
red subpopulations respectively, but only 61 of 77 GCs for 
the green subpopulation as we restrict the green sample to 
include only GCs with galactocentric radii smaller than 5 
arcmin. We impose this galactocentric radius restriction on 
the green subpopulation because beyond this radius it con- 
tributes less than ten per cent of the total GC population. 
The GC subpopulation kinematics are fitted as a function 
of galactocentric radius in rolling bins of 22 GCs and the re- 
sults of the kinematic fits are shown in Fig. [9] The bin size of 
22 GCs was chosen to maximize the radial range of the kine- 
matic analysis. This is significantly smaller than the bin size 
used for Fig. [7| because inter subpopulation contamination 
has been reduced. We fit the kinematics in the case where 
the PAfci„ is allowed to vary with galactocentric radius. 

The blue subpopulation shows significant rotation at 



all radii. There is an indication that the PAtin of rotation 
twists from ~ 200° to ~ 310° between 3 and 7 arcmin (20 
and 47 kpc). The rotation velocity is consistent with being 
fiat and Vrot ~ 100 km s^^ over the entire radial range of 
the blue subpopulation. The velocity dispersion of the blue 
subpopulation is also consistent with being fiat and a ~ 200 
km s~^. The radial kinematic profiles are consistent with the 
values found for the blue GCs as seen in Fig. [7] as well as the 
fitted values for the entire blue subpopulation sample (Ta- 
ble [l]). The kinematics of the blue GC system varies more 
significantly with galactocentric radius than GC colour and 
by removing most of the contaminating GCs coming from 
the green subpopulation as well as dividing the blue GC by 
radius, we achieve tighter constraints on the subpopulation 
kinematics. The twist in the PAkin could suggest that the 
radial outskirts of the blue subpopulation have been affected 
by accretion of galaxies and associated GCs, or a tidal in- 
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Figure 10. Kinematic parameters of GCs as a function of galac- 
tocentric radius. Rotation parameter Vrot/o" (top panel) as well 
as the root mean squared velocity Vrms (bottom panel) for GC 
subpopulations potted against galactocentric radius. Linestyles 
are the same as those in Fig. [9] This figure shows the case where 
the position angle is allowed to vary. The rotation bias is small 
everywhere except inside a galactocentric radius of 2 arcmin for 
the red subpopulation. Here the plotted value for VVot/f is sig- 
nificantly overestimated but because there is no fitting involved 
in obtaining the Vrms this value is not affected. 



teraction. This kinematic twist could also be a signature of 
triaxiality in the blue GC subpopulation. 

The green subpopulation rotates relatively fast at small 
radii {Vrot ~ 200 km s~^) but the rotation decreases quickly 
with increasing galactocentric radius to ~ 75 km at 3.5 
arcmin. The radially fitted PAfei„ is consistent with being 
flat (150 ± 30°) and consistent with the PAfci„ fitted to the 
entire subpopulation (144±23°). The velocity dispersion of 
the green subpopulation also decreases with galactocentric 
radius from 260 to 180 km s~^. 

The red subpopulation rotates fastest between 2 and 4 
arcmin, with VKkin ~ 222° . The rotation velocity amplitude 
here is Vrot ~ 180 km s~^.The fitted PAfci„ within 2 arcmin 
is not as well constrained as that fitted between 2 and 4 
arcmin, and within 2 arcmin the rotation parameter fails 
the criterion defined in Equation [5] The rotation of the red 
GC subpopulation at small galactocentric radii is consistent 
with zero rotation and the value of ~ 75 km s~^ fitted here 
is due to rotation bias caused by the free variation of PAfci„. 
The velocity dispersion of the red GCs decreases with radius 
from a ~ 220 to ~ 150 km s~^. The red subpopulation 
is dispersion dominated in the central regions and rotation 
dominated in the outer regions of the galaxy. 

Plotting the rotation parameter CK ot/"") in Fig- 



10 



see that the red subpopulation shows very different proper- 
ties to the blue and green subpopulations. This ratio is one 
measure of the diskyness of the system, where Vrot /o > 1 
can indicate a disky system. It is only the red GC subpop- 
ulation that shows rotationally dominated characteristics, 
and only at galactocentric radii beyond 2.5 arcmin, where it 



rotates with a PAfei„ of 
NGC 4365). 

The root 



222° (photometric major axis of 

nean squared velocity 

shows variation with radius 



Vrms ' 

for all three subpopulations but changes most significantly 
in the green subpopulation. The green GC subpopulation 
Vrms decreases from ~ 270 to ~ 190 km s~^ over 1.5 
arcmin in galactocentric radius. The sharply declining 
Vrms seen in the case of the green subpopulation could 
be explained by a truncation in their density distribution 
or a bias towards radial orbits relative to the other GC 
subpopulations. 



4 DISCUSSION 

Since most giant elliptical galaxies contain only two metal- 
licity subpopulations of GCs it is likely that one of NGC 
4365's three GC subpopulations is the result of an addi- 
tional and unusual formation/accretion process. Determin- 
ing which of the GC subpopulations is additional may indi- 
cate which mechanism caused it and could be a key to un- 
covering the formation mechanisms of GC subpopulations 
in ordinary elliptical galaxies. We summarise the kinematic 
and photometric properties of the three GC subpopulations 
as well as the two main components of the galaxy starlight 
(the KDC and the bulk of the stars) in Table|2]and use these 
properties to discuss possible formation scenarios. 

We find three distinct kinematic signatures for the three 
GC subpopulations. The blue subpopulation rotates at all 
observable radii and its kinematic position angle twists 
slowly with increasing galactocentric radius from 200° to 
312°. Its rotation is unrelated to NGC 4365's stellar light. 
The green subpopulation rotates with the bulk of NGC 
4365's stellar fight along the photometric minor axis of the 
galaxy (132°). The red subpopulation rotates only at large 
radii along the photometric major axis (222°) of NGC 4365 
but in the opposite direction to the KDC of the galaxy (38°). 
The distinct kinematics of the red and green GC subpopula- 
tions, seen clearly in Fig. |ll[ provide evidence that the green 
GCs form a separate subpopulation and are not simply the 
metal poor tail of the red subpopulation. 



The Blom et al. (2012 I photometric analysis found that 
the radial density distribution and colour of the red GC sub- 
population most closely matched the stellar light of NGC 
4365 and, while the axis ratio of the green subpopulation 
more closely matched that of the starlight, they concluded 
that there was evidence that the red subpopulation was re- 
lated to the bulk of the galaxy starlight. The analysis here 
shows that the rotation position angle, and possibly also the 
rotation velocity amplitude, of the green GC subpopulation 
most closely match the extrapolated stellar kinematics. Thus 
it is not immediately clear which of the two is the additional 
subpopulation. 

If the red GC subpopulation is additional it could have 
been created by an additional, later (still > 10 Gyr ago) 
merger ( |Zepf fc Ashman|1993[ ). This might explain why the 
red subpopulation rotation is decoupled from the bulk of the 
starlight and the green subpopulation rotates with the bulk 
of the stars. The problem arises when comparing this red 
subpopulation to those of other galaxies of similar luminos- 
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I 2012 I, * from Krajnovic 



Table 2. Su mmary of the NGC 4365 galaxy system pro perties. Values mark ed with t are taken from 

|et al.| | |201l] l, * from |Ferrarese et al.| l |2006[ l and * from |Davies et al.| pOOlj l, while all other values are derived in this work. Notes: ^The 
position angles for the blue and green subpopulations were qualitatively similar to the galaxy light and were therefore assumed equal to 
it. ^This effective radius measurement is not obtained directly from a Sersic profile fit to the green subpopulation but by comparing the 
power law slope of the green and red subpopulations. ^The g' — z' value was converted to a g' — i' colour value. ''These are maximum 
velocities observed in a limited field-of-view not mean or median velocities. ^This is the largest radial extent of the KDC and not its 
effective radius. 





g' - i' (mag) 


"^^photom ( ) 




Vrot (km s 1) 


(7 (km s ^ ) 


Qphot 


_Re (arcmin) 


Blue GCs 
Green GCs 
Red GCs 

Stars 

KDC 


0.80 t 
0.97 t 
1.13 t 
1.15 t 
1.40 


~ 42/222 t.i 
~ 42/222 t.i 

undefined 
42/222 ± 3 t 
40.9 ± 2.1 « 


269±2| 
144±ii 
194±23 
145 ±6.5 " 
38±5 • 


67±t| 
81±39 
94±|| 
60.9 "'^ 
63 ±8 


178±10 
217±io 
202±15 
260 ± 10 
262 ± 3 


0.56 ±0.08 t 
0.70 ± 0.08 t 
0.97 ± 0.08 t 
0.75 ±0.03 t 
0.75 ±0.08 t 


7.30 ±0.68 t 

~ 2.97 t>2 
3.17± 0.19 t 
2.1 ± 0.11 t 
0.12 ±0.02 ''S 



ity to NGC 4365, as it is not an outlier but consistent with 
the colour of red subpopulations of other galaxies of similar 
luminosity/mass ( Peng et al.|2006 l. 

If the green GC subpopulation has been recently do- 
nated from a medium sized galaxy that accreted on to NGC 
4365 then the remnant of this accreted galaxy might still be 
visible in the form of a stellar stream across NGC 4365 as re- 



ported by |Bogdan et al.| ( |2012[ ). |Blom et al.| ( |2012[ ) reported 
that the green subpopulation is centred aX g' — i ~ 0.97 and 
we can estimate the absolute magnitude range for the ac- 
creted galaxy by comparing the peak colour of the green sub- 
population with the relationship between peak GC subpop- 
ulation colour and galaxy luminosity empirically determined 



by Peng et al. (20061. In this case the green subpopulation 



of NGC 4365 is the original metal-rich (red) GC subpopu- 
lation of a —14 < Ms < —15 dwarf galaxy and we expect 
to find primarily blue and green GCs in the stellar stream 
of NGC 4365. The total stellar light observed in the low 
surface brightness stream is _B ~ 13.5 mag (Bogdan et al 



2012 this a minimum brightness for the accreted galaxy) 
corresponding to Ms ~ -18.3 at the distance of NGC 4365. 
There is a sizeable mismatch between the predicted abso- 
lute magnitude for the dwarf galaxy and the much brighter 
integrated stellar light from the stream, suggesting that the 
stellar stream is not the remnant of the galaxy associated 
with accreted green GCs. This explanation has several other 
problems. It is still to be explained why the accreted green 
subpopulation is so centrally concentrated while the stellar 
light of the accreted galaxy is still on the outskirts of NGC 
4365. Also why does the green subpopulation rotate with 
the bulk of the stellar light of NGC 4365 (along the minor 
axis of the galaxy) rather than along the stream orientation, 
which is roughly aligned with the major axis of NGC 4365? 

If the accretion event which brought the green GCs into 
NGC 4365 was not recent, the currently observable stellar 
stream may be unrelated to the green GCs. We can still 
investigate this claim based on the number of green GCs (see 
a similar analysis by Romanowsky et al.|2012 l. Blom et al 



( 2012 1 reported that the green subpopulation contributes 17 



percent of the total GC population or 1100±20 GCs to NGC 
4365's GC system (Nrot = 6450±110). We can calculate the 
absolute magnitude range of the accreted galaxy using the 
observed relations hip between GC specific frequ ency (Sjv = 
A^GclO° '*'*^^+^''MHarris & van den Bergh|l98lk and galaxy 



luminosity ( Peng et al.|2008 1 . To be consistent with the S n 
range observed for GC systems in the literature, the accreted 
galaxy must be brighter than Mv ~ —20.5 [Mb ~ —19.5). 
The number of GCs in the green subpopulation far exceeds 
the expected number of GCs that would be associated with a 
dwarf galaxy. If the green subpopulation was donated by an 
accreted galaxy (at any stage in the history of NGC 4365) 
the mean colour of the green GCs suggests a very small 
galaxy whereas the number of green GCs suggests a large 
galaxy. 

If the green GCs were accreted from a group of dwarf 
galaxies it might account for the large number of GCs with 
colours that indicate a small host galaxy, but the group of 
dwarfs would need to contain many tens of members to ac- 
count for the large number of inferred green GCs. 

We also consider the possibility that both the red and 
green subpopulations were formed in-situ. In this case the 
green subpopulation would have formed unusually early in 
the dense core of the galaxy while the red GCs formed later 
when the core of the galaxy was slightly bigger and slightly 
more metal rich. This scenario might explain why the green 
subpopulation is more centrally concentrated than the red 
subpopulation and why the green GCs rotate with a signifi- 
cant fraction of the starlight. It does not explain why we do 
not detect a fraction of the stars rotating with the red GCs 
or what would have caused two epochs of early in-situ GC 
formation. 



Bogdan et al. ( 2012 1 s tudied the lenticula r galaxy NGC 
Grogin et al.||l998 l embedded in 



4342 (V;^^ = 751 km s" 

the stellar stream across NGC 4365 (see Fig. |ll[ ) and inferred 
a massive and extended dark matter halo around NGC 4342. 
The central velocity dispersions of NGC 4365 (261 ± 7 km 
s"^) a nd NGC 4342 (244 ± 1 1 km s"^) are consistent within 
errors ( Bernardi et al.|2002 I implying similar mass galaxies. 



It is possible that this very massive, if relatively faint, galaxy 
could have passed close by NGC 4365 and dragged some 
of NGC 4365's own stars out from the galaxy to form the 
stream. It is also possible that the stream stars have been 
stripped from NGC 4342 itself. These scenarios both agree 
with our observation of low velocity GCs (at ~ 700 km s~^, 
similar to the Vsys of NGC 4342) in NGC 4365 (Fig. |4| that 
are aligned with the stellar stream (Figs[6]and |ll[ ). The low 
velocity GCs are consistent with the colour distribution of 
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Figure 11. Schematic representation of the NGC 4365 stellar 
light and GC system. The diagram is oriented with North up and 
East left. Positions of low velocity GC outliers are shown with 
black dots. Ellipses show the 2Re extent of the galaxy light (black) 
and three GC subpopulations in blue, green and red respectively. 
The shape and relative size of the stellar light stream (reported 
in |Bogdan et al.|2012[ | is shown with a grey irregular outline. The 
rotation direction of the galaxy light as well, as the green and red 
GC subpopulations, are shown in the top right of the diagram. 
The rotation direction of the blue GCs is not shown because there 
is a significant twist with galactocentric radius. A cross marks the 
position of NGC 4342. The low velocity GCs are generally aligned 
with the stellar stream and have a similar velocity to NGC 4342. 



NGC 4365 's underlying GC population (Fig.[5| but it is still 
unknown how their colours compare to those of NGC 4342. 

It is tempting to suggest that the close passage of NGC 
4342 might have stretched NGC 4365 starlight and GC sys- 
tem into its current shape. If NGC 4365 was round or only 
mildly elliptical (with a PAphot close to the current PAkin 
of the bulk of stellar hght) before NGC 4342 made its close 
pass we would expect to see remnants of the original shape in 



the central regions of NGC 4365. |Ferrarese et al. ( 2006 1 find 
the galaxy ellipticity drops below 0.05 in the g and z bands 
while the PAphot twists to ~ 90° at radii smaller than 0.3 
arcsec. These photometric properties only appear at scales 
smaller than the KDC of NGC 4365 and are unlikely to re- 
late to the bulk of NGC 4365. This scenario might explain 
the anomalous 'rolling' starlight as an artefact of very recent 
perturbations to NGC 4365's photometric shape but would 
not explain the presence of NGC 4365's KDC . The KDC 
could be a remnant of an old merger, as its age is observed 
to be > 10 Gyr and consistent with the galaxy starlight 



(Davies et al. 20011 or the consequence of our particular 



van den 



viewing angle into an intrinsically triaxial system 
Bergh et al.|1991 |. 

Fig. |11| summarises the photometric and kinematic in- 
formation of the NGC 4365 system graphically. The blue 
GC subpopulation is the most radially extended GC sub- 
population and is more elongated than the galaxy light. Its 
rotation shows a slow twist from alignment with the galaxy 
major axis in the inner parts to the galaxy minor axis in the 



outer parts. The amplitude of its rotation and velocity dis- 
persion are constant with radius. The green subpopulation is 
the most centrally concentrated, roughly as elongated as the 
galaxy light and shows 'rolling' kinematics like the galaxy 
light {PAkin ~ galaxy minor axis). Its rotation and velocity 
dispersion amplitude decrease with increasing radius. The 
red subpopulation is almost circular and rotates only in the 
outer parts, along the galaxy major axis (in the opposite di- 
rection to the KDC). The low velocity GCs are aligned with 
the stellar stream and contain GCs of all colour subpopu- 
lations. The kinematics of the galaxy and GC system are 
complex but it appears that the stellar stream is not con- 
nected to the presence of any one of the GC subpopulations. 



5 CONCLUSIONS 

NGC 4365 is a giant elliptical galaxy with rare stellar kine- 
matic properties and also, very unusually, three globular 
cluster (GC) subpopulations. Recently, a stellar stream run- 
ning across NGC 4365 and its nearby neighbour, NGC 4342, 
was discovered. 

With high velocity resolution Keck/DEIMOS spectra 
we analyse the kinematics of a large number and spatially 
extended sample of GCs around NGC 4365. These data al- 
low us to separate the GC system into its three colour sub- 
populations and analyse the kinematics of each subpopula- 
tion in detail. We find distinct rotation properties for each 
GC subpopulation. This indicates that the three GC sub- 
populations are distinct and that the formation history of 
NGC 4365 is complex. We conclude that it is unlikely that 
the third (green) GC subpopulation is related to the exis- 
tence of the stellar stream and that it is also unlikely to have 
originated from an accreted galaxy. 

We also find a further group of low velocity GCs (cov- 
ering the full range in colour), which might be related to 
the stellar stream extending across NGC 4365 towards NGC 
4342 (~ 35 arcmin away at a similar systemic velocity to 
these GCs). Future analysis of spectroscopic observations of 
GCs in the stream will further illuminate the complex for- 
mation history of NGC 4365. 
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APPENDIX A: TOTAL SUBPOPULATION 
KINEMATIC FITS 

The colour peaks and widths of the blue, green and red sub- 
populations cause significant overlap in the colour ranges of 
the subpopulations. We explore the effect of splitting the 
populations in different ways to determine the best way to 
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Table Al. Results for kinematic fits to various divisions of the three GC subpopulations. We tabulate the rotation velocity {Vrot), 
velocity dispersion (tr) and kinematic position angle (PAj.j„) for each subpopulation in all the split definitions. Results are shown where 
the axis ratio {quin = Qphot) is fixed to the photometric values for each subpopulation (qbiue = 0.56, qgreen = 0.70, qred = 0.97 and 
qgal = 0.75) and where it is fixed to 1. Numbers of GCs in each subpopulation representative sample are given in brackets after the 
definition parameters. The major axis of the NGC 4365 stellar light is 42° or 222° . Probability is denoted with Pcolour 



Split 
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Defn 

of qkin 


Vrot 

(km s-1) 
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minimize the contamination between samples of the subpop- 
ulations and maximize the number of GCs in each subpop- 
ulation sample. 

The first subpopulation split definition is based on the 
probability that each GC is assigned of being in either the 



blue, green or red subpopulations by Blom et al. (20121 



They assigned each GC a probability of belonging to the 
blue, green and red subpopulations {Pbiu(; + Poreen + Pred = 
1) based on its colour and galactocentric radius. We define 
medium probability samples with P > 0.80 to represent the 
three subpopulations and find 53 blue, 53 green and 65 red 
GCs with line-of-sight radial velocity measurements. These 
samples are mostly free of contamination from the other 
subpopulations but have a relatively small number of GCs 
to be analysed for each subpopulation. 

The next subpopulation split definitions cut the GCs 
into three subsamples at different g' — i' colours. For the first 
colour cut split definition we split the sample at g' — i' = 0.9 
and g' — i' = 1.1 where we see dips in number of the photo- 
metric sample colour distribution (see Fig.[2|. This subpop- 
ulation split definition divides the kinematic sample roughly 
equally between the three subpopulations. The next two 
colour cut split definitions widen the green sample defini- 
tion first to redder colours (colour cut 2) and secondly to 
bluer colours (colour cut 3). The repeated fits for the blue 
(in colour cut 1 and colour cut 2) and red samples (in colour 
cut 1 and colour cut 3) show the level to which results are 
affected by the randomisation in the bootstrap technique of 
determining uncertainties. 

We also show the kinematic fits for both cases of the 
kinematic axis ratio. In the first case qtin is fixed to the pho- 
tometric value for each GC subpopulation (qkin = Qphotom 
from Blom et al. 2012[ ) and in the second the axis ratio is 



fixed to unity [qkin = !)• We do not find any significant dif- 
ferences in the fitted kinematics between the two axis ratio 
cases and conclude that this kinematic fitting method is not 
very sensitive to the kinematic axis ratio. In all split defini- 
tions for all subpopulations, the kinematic fits for the two 
cases are consistent within la uncertainty and the computed 
A/ndf values are equal. All further kinematic fits are done 
in the case where qkin = 1 as this is the more general case. 

The best kinematic fits are found for the blue GCs in 
colour cut 3, green GCs in colour cut 1 and red GCs in colour 
cut 2. We determine this best kinematic fit to be the subpop- 
ulation sample has the most constrained position angle (for 
q=l), see Table Al and lowest value of A/ndf, see Table 



|A2| The blue subpopulation representative sample in colour 
cut 3 A/ndf is marginally larger than in colour cut 1 or 2 but 
the position angle range is significantly smaller (62°) than in 
either colour cut 1 or 2 (90° and 98° respectively). The green 
subpopulation representative sample position angle range in 
colour cut 1 (50°) is only slightly smaller than that of colour 
cut 2 (52°) but the A/ndf is significantly smaller. The red 
subpopulation representative sample A/ndf is the same for 
each colour cut but the position angle range is significantly 
smaller in colour cut 2 (47°) than either of the other colour 
cuts (91° and 101°) where the PAkin is not constrained to 
< 90°. 



APPENDIX B: RADIAL KINEMATIC FITS 
WITH FIXED POSITION ANGLE 

Here we compare the radial kinematic fits for the three sub- 
populations in the case where the PA is fixed to the value fit- 



ted for the whole subpopulation (Fig. Bl I to the case where 
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Table A2. Minimisation parameter divided by degrees of freedom 
(A/ndf) for kinematic fits to various divisions of the three GC 
subpopulations. Values for the best colour cuts are highlighted. 
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Figure Bl. Kinematics as a function of radius for NGC 4365's 
three GC subpopulations. GC subpopulations are defined using 
the best colour cuts highlighted in Table |A1| and the PA^^^ is 
fixed to the values found there. Symbols are as in Fig. [9] 



the PA is allowed to vary with radius (Fig. [9|. The fitted 
values for a do not change at all when the position angle of 
rotation is fixed. The uncertainties in Vrot, are larger when 
PAfcin is fixed than when the PAfci„ is allowed to vary. If the 
reduction in number of free parameters, in the case where 
PAkin is fixed, does allow a more stable fit to the Vrot and 
a we would expect the uncertainties (68 per cent confidence 
intervals) to be smaller there. Since the opposite is true it is 
more likely that the fit is better where the PAkin is allowed 
to vary, and since the PAkm is well constrained at almost all 
radii in that case it is also likely that the Vrot is also more 
reliable than the case where PAfei„ is fixed. 

The one exception is for red GCs less than 2 arcmin 
from the galactic centre. The varying PAfei„ is not well con- 
strained there and rotation parameter (K-ot/o") fails the cri- 
terion in Equation [s] ( [Strader et al.|2011 1. In the case where 
PAkin is fixed Vrot is correctly fitted to be consistent with 
zero. The decrease in Vrot at certain other points in Fig. |B1| 
(e.g. ~ 3 arcmin for the blue subpopulation) is due to the 
offset between the value the PAkin was set at and the PAkin 
at which the system is actually rotating at that radius. 
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